Diffusion Limited Aggregation with Power-Law Pinning 
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Using stochastic conformal mapping techniques we study the patterns emerging from Laplacian 
growth with a power-law decaying threshold for growth ii^^ (where -Rjv is the radius of the A''— 
particle cluster). For 7 > 1 the growth pattern is in the same universality class as diffusion limited 
aggregation (DLA) growth, while for 7 < 1 the resulting patterns have a lower fractal dimension 
,0(7) than a DLA cluster due to the enhancement of growth at the hot tips of the developing pattern. 
Our results indicate that a pinning transition occurs at 7 = 1/2, significantly smaller than might 
be expected from the lower bound Qmin — 0.67 of multifractal spectrum of DLA. This limiting 
case shows that the most singular tips in the pruned cluster now correspond to those expected for 
a purely one-dimensional line. Using multifractal analysis, analytic expressions are established for 
D{"f) both close to the breakdown of DLA universality class, i.e., 7 < 1, and close to the pinning 
transition, i.e., 7 > 1/2. 



B 

o 
o 



> 
00 
00 

o 
m 
o 



X3 
O 

o 



X 



PACS numbers: 05.45.Df, 61.43.H\ 



I. INTRODUCTION 



As we will show below, as 7 decreases from 7 > 1 to- 
ward a critical value 7 = 1/2 (which corresponds to the 
most singular possible behavior for the Laplacian field, 
that occuring at the tip of a line), there is a continuous 
transition from DLA toward lower dimensional shapes for 
which the multifractal spectrum is necessarily different 
from that of DLA. We study this transition in terms of 
the fractal dimension D('y) of the emerging patterns, and 
we derive analytic expressions for the behavior of £'(7) in 
the range 7^1 where the DLA universality class breaks 
down and 7 > 1/2 close to the pinning transition. 



Nonequilibrium growth models leading naturally to 
self-organized fractal structures, such as diffusion lim- 
ited aggregation (DLA) ^ , are of continuing interest due 
to their relevance for many important physical processes 
including dielectric breakdown Q, electrochemical depo- 
sition U3 , and Laplacian flow [5| . 

A powerful method, namely iterated stochastic confor- 
mal mapping 0,0) has been already successfully appHed 
to generate and analyze DLA 's', "yl and Laplacian 
growth patterns in two dimensions. This has provided an 
alternative way to address many of the important open 
questions related to pattern formation in DLA in two 
dimensions, one of these being the existence of minimal 
fields for growth at the boundary of the growing cluster. 
In previous work [Tll | we studied the properties induced 
by a fixed, material dependent critical field Ec for growth, 
and showed that in the presence of such a threshold all 
clusters ultimately become pinned and, in addition, this 
simple constraint has remarkable consequences for the re- 
sulting patterns - the rich, branched structure of DLA is 
replaced by a much lower dimensional shape consisting 
of a few surviving branches. 

In this paper we address a similar, but significantly 
more important question because of its relationship to 
the multifractal spectrum of DLA, that of what happens 
when there exists a critical field for growth on the bound- 
ary of the cluster, field decaying like RJp as the cluster 
increases in size. We shall call this model the "7-model". is the radius of the N particle cluster, i.e.. 



II. MODEL AND THEORETICAL 
BACKGROUND 

Witten and Sander have shown that the growth 
probability at any point s on the boundary of a DLA 
cluster of length L is given by the harmonic measure 
P(s) = |(VF)(s)|/ ds'\{VV){s')\, where V{r) obeys 
Laplace's equation V^F — subject to the boundary 
conditions V = on the (evolving) boundary of the clus- 
ter and y '-^ In r as r — > cx) (corresponding to a uniform 
flux of particles far away from the cluster). 

The model we study is a variant of the two-dimensional 
DLA growth model described above in which growth is 
disallowed at points on the cluster boundary where the 
probability for growth is smaller than a critical value 
^N^i where Rn (the exact meaning will be defined later) 
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where a(s') is the multifractal exponent at point s' on the 
cluster boundary, L is the length of the boundary, and the 
step function 9{-f — a{s')) ensures that only those regions 
of the cluster boundary obeying |Vy(s')| > RJf con- 
tribute to the normalization integral. Estimates of this 
integral will be very important in our analysis of the frac- 
tal dimension of the growing "7-cluster" . Since we need 
to calculate the harmonic measure on a freely evolving 
interface, this is handled by using conformal mapping 
techniques 0, 13. The method was presented in great 
detail in Refs. [a, Q, and thus here we will just briefly 
review the main results. 

The basic idea is to follow the evolution of the con- 
formal mapping z = (f>^"-'(aj) of the exterior of the unit 
circle in a mathematical tj-plane onto the complement of 
the cluster of n particles in the physical z-plane rather 
than directly the evolution of the cluster's boundary. The 
initial condition is chosen to be <i>^°''(a;) = ui. The pro- 
cess of adding a new "particle" of constant shape and 
linear scale ^/Xo to the cluster of (n — 1) "particles" at 
a position s which is chosen randomly according to the 
harmonic measure is described via a function (pxfii'^), 
where 
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which conformally maps the unit circle to the unit circle 
with a burnp of linear size \/X localized at the angular 
position 6 lf||. The shape of the bump depends on the 
parameter a. Following the analysis in we have used 
a — 0.66 througout this paper, as we beheve the large 
scale asymptotic properties will not be affected by the 
microscopic shape of the added bump. 

The conformal map for an n-particle cluster 
can be built by adding one "particle" to an (n— l)-particle 
cluster $'^"^-'^^(0;), resulting in the recursive dynamics 



(3) 



which can be solved in terms of iterations of the elemen- 
tary bump map 4>x^ q^{ijj), 

= 0Ai,0i o o • ■ • o 0A„,e„(t^) • (4) 

In Eqs. El and 01 the angle 0„ G (0, 27r] at step n is ran- 
domly chosen since the harmonic measure on the real 
cluster translates to a uniform measure on the unit circle 
in the mathematical plane. 
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is chosen in order to ensure that the size of the bump in 
the physical z plane is a/Aq- Since ^/Xq is a natural length 
scale in the problem, in can be scaled out by measuring 
all the lengths in terms of it. We re-emphasize here that 
although the composition Eq. 01 appears at first sight to 
be a standard iteration of stochastic maps, this is not so 
because the order of iterations is inverted - the last point 
of the trajectory is the inner argument in this iteration. 
As a result the transition from $^"~^)(a;) to is 
achieved by composing the n former maps Eq. ^starting 
from a different seed. Finally, identifying 8] the radius 

(n) 

Rn of the growing pattern with the coefhcient F{ = 
n"=i(l + Xi)°- in the Laurent expansion of 
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the constraint to grow only at values of 
Eq. m translates into 



1 



■ , (7) 
which obey 

(8) 



This constraint is implemented as follows. At step n. On 
is choosen from a uniform distribution in (0,27r], inde- 
pendent of previous history. If it obeys the constraint 
given by Eq. [HI accept this value of 0„, otherwise repeat 
until the constraint is obeyed. 



III. RESULTS AND DISCUSSION 

For the model defined in Sec. II one would expect the 
resultant patterns to have fractal shapes which depend on 
7, and in order to characterize these shapes we will focus 
on the scaling behavior of the first Laurent coefficient 
f["'\ Following the arguments in Ref. Q, for a given 
value 7 one expects a scaling law of the form 



(n) 



,1/-D(7) 
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where -D(7) is the effective fractal dimension of the re- 
sulting cluster. 

We have simulated the model defined in Sec. II for 
a number of values 7 in the range 1/2 < 7 < 1.2 and 

we have calculated F^^"^ as an average over 100 clusters 
(for 7 > 0.65), and, respectively, over 20 clusters (for 
7 < 0.65) of size N = 40000. In Fig.Uwe show typical 
clusters of size N = 40000 for 7 = 1.20, 0.75 and 0.55, re- 
spectively. It can be easily seen that the rich, branched 
structure of the cluster at 7 = 1.20 changes, as 7 de- 
creases toward 7 = 0.55, into a much lower-dimensional 
shape with only a few branches surviving. This can be 
intuitively understood by considering the effect the pin- 
ning probability R~'^ has on the multifractal spectrum 
of the cluster. From multifractal scaling \12L |lj| we 
know that the interface of a fully developed DLA cluster 
consists of sets of Ndla^cx) ~ i?;^"^^^"^ sites with grow- 
ing probabilities \W\ ~ and we will assume that 
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FIG. 1: Typical clusters (size = 40000) grown with (a) 7 ■ 
1.20, (b) 7 = 0.75, and (c) 7 = 0.55, respectively. 



such a structure is also valid for the clusters grown with 
pinning probability R^"'- As growth proceeds, the lowest 
probability sites (large a), normally deep in the fjords, 
will be pinned first, with the hot tips surviving longest, 
leading to "pruning" of the branches where the tips have 
a singularity a > 7. 

As anticipated, for all the values 0.5 < 7 < 1.2 that we 
have tested the coefficient F^"' has a clear power-law de- 
pendence on the size n, as shown in Fig.|2a). Assuming 
the exponent 1/D{'j) to be related to a fractal dimen- 
sion as given by Eq. |3 the dependency I? (7) (shown in 
Fig.|2Ib)) is obtained from a power-law fit to the data. It 
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FIG. 2: (a) F^'"' as a function of n for clusters grown with 
7 = 1.20, 0.75, 0.65, and 0.55, respectively (log-log plot). Also 
shown (dashed lines) are the limit cases of a DLA cluster, 
F^"'' ~ n^/-^"^^, where Ddla = 1.71, and of a line cluster, 
_F^'"' ^ n, respectively, (b) The effective fractal dimension 
-D(7) obtained from F^' - ni/°(^) as a function of 7 (cir- 
cles). The dotted line is just a guide to the eye. The solid line 
is the theoretical prediction (Eg. I12II with f-,{'y) — foLA^'y) 
from Ref El. 



can be seen that at values of 7 > 1 the behavior is close 
to that of a DLA cluster, i.e., £'(7 > 1) — > Ddla — 1.71, 
while for 7 > 1/2 the behavior of the radius f["'^ tends to 
n, i.e., D(7 1/2) 1, thus the behavior of a growing 
line. 

In order to understand these results theoretically let us 
begin with a very simple argument based on the assump- 
tion that the clusters have a multifractal spectrum in the 
sense described above, in that the interface consists of 
sets of A/'-y(a) ~ Rn^°'^ sites with growth probabilities 
|Vy| ^ i?,7". Since the constraint will cut-off growth at 
regions in the cluster with exponents in the range a > 7 
of the multifractal spectrum, we can write down the fol- 
lowing equation for the rate of growth of the cluster in 
the presence of the barrier in terms of the rate of growth 
of a DLA (no barrier for growth) cluster 
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where f{amin) = 0. The enhancement of growth comes 
from Eq.^ together with the estimate 

(10) 

Now, we know that the harmonic measure is concentrated 
at a = 1, and thus for 7 > 1 the integral is dominated 
by the value of the integrand at a = 1 , while for 7 ^ 1 it 
is dominated by the value at 7, and thus 



dR 

In 



for 7 > 1, 



X < ' "t , \ ' ' (11) 

DLA [R-'-f-^-'K for7<l, 



where if (7) is some constant independent of N. Eq. [TTl 
therefore implies 



A 

DLA 



+ 77(7) -7, 



for 7 > 1, 
for 7 < 1. 



(12) 



For the case where 7^1, i.e., close to the breakdown in 
the DLA universality class, we may assume that the mul- 
tifractal spectrum of the cluster is only weakly perturbed 
from its value in DLA and therefore f-yi'^) ~ foLAil), 
or 



[DdLA + JDLA[1) 



for 7 > 1, 
for 7 < 1. 



(13) 



This prediction may be tested using the recently com- 
puted fDLA{ot) spectrum 'l£?|. The results are shown in 
Fig-EIb) (solid line), and it can be seen that for 7 ~ 1 
the theoretical predictions are indeed close to the mea- 
sured values 13(7). The discrepancies at smaller values 
of 7 can be attributed to the fact that the multifractal 
spectrum of the cluster is not exactly the one of DLA, 
and when 7 1/2, as we will now discuss, it actually 
may be expected to deviate significantly from that of a 
DLA. 

For 7 > 1/2, the change in the multifractal spectrum 
from that for DLA is significant. For example, it is known 
[l5lll6^ that for the DLA spectrum amin — 0.67, while we 
see that growth continues significantly below this value, 
with amin = 1/2 being the asymptotic limit. Our simu- 
lations show strong evidence for this limit, as can be seen 
in Fig.|2Ia): the average number of attempts for growing 
clusters of sizes N = 20000 and N = 40000 (scaled by 
the actual size of the cluster) exibits a steep increase as 
7 0.5. 

We shall assume that for these highly pruned "7- 
clusters" there is a well defined limiting form for the 
multifractal spectrum, f-^f{a), defined for 1/2 < a < 7 
and obeying fj{l/2) = 0, which is a monotonically in- 
creasing function of 7. Defining 7(7) = f-y{l), we can 
use a reasoning similar to that used by Turkevich and 
Sher in their estimate of the fractal dimension of 
DLA to estimate /(7). The idea is that for 7 > 1/2 only 
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FIG. 3: (a) Average number of attempts for growing clus- 
ters of sizes iV = 20000 (dashed line) and iV = 40000 (solid 
line), scaled by the actual size of the cluster, respectively, as 
function of 7 (b) The multifractal spectrum f('y) from Eg. 1151 
(symbols) in comparison with the multifractal spectrum of 
DLA (solid line) from Ref. E^. 



the "hottest" tips contribute to growth, and thus one can 
write 



dR 
dN 



^-l/2-/(7)+7^ 



(14) 



where the last relation follows from the fact that the in- 
tegral is dominated by the value of the integrand at 7. 
Thus, we obtain 



/(7) = i?(7)+7-3/2. 



(15) 



Since from simulations we know the values D{^j), Eg. 1151 
(which is valid for 7 > 0.5) allows the calculation of 
the upper multifractal exponent /(7). As shown in 
Fig. Olb), the pinning threshold R '' leads to a shift- 
ing of the multifractal spectrum for a < 7 to the left, 
i.e., f-/{a) > foLAict) (more hot tips, and larger fields at 
those hot tips, due to pruning). 

These results can be intuitively understood as a flow of 
singularities away from 7 (which acts as an unstable fixed 
point of the dynamics). For any particular value ao < 7, 
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what happens while the cluster evolves is that screening 
is reduced compared to DLA and therefore there is a flow 
of singularities ao — > ai with ai < aQ. In addition, new 
singularities with a < cx'min^K can be created. Thus, we 
would expect that the number of singularities A/'-y(ai) « 
J^DLAiao) or 

/^(ai(Q!o)) w /uLA(ao)- (16) 

On the other hand for ag > 7 the singularity flow 
cko ^ cti can only act toward an increase ai > ao since 
such points can never grow and thus can only either keep 
their original singularity or get a higher value of a during 
growth. 

IV. CONCLUSIONS 

Using the stochastic conformal mapping techniques 
we have studied the patterns emerging from Laplacian 
growth with a power-law decaying threshold for growth 
^jv^: 7 ^ 1/2. We have shown that due to the enhance- 
ment of growth at the hot tips as 7 decreases the growth 
evolves from patterns in the DLA universality class for 
7 > 1 to clusters with a lower fractal dimension I? (7) for 
7 < 1 due to the enhancement of growth at the hot tips. 
We have presented evidence that 7 = 1/2, corresponding 
to the singularities at the tip of a purely one-dimensional 



(line) growth pattern, is the lower limit for growth, with 
all clusters becoming ultimately pinned for 7 < 1/2. By 
using multifractal analysis, we have proposed analytic ex- 
pressions for -D(7) for both 7^1 near the breakdown of 
the DLA universality class and near the pinning transi- 
tion 7 > 1/2. Finally, we have shown that in the small 
7 range the multifractal spectrum of the resulting clus- 
ter is significantly changed from that of a DLA. We have 
suggested that this change may be due to a flow of sin- 
gularities with 7 acting as an unstable fixed point of the 
dynamics, but further work will be necessary to fully elu- 
cidate this point. 
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